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Abstract 

Three methods: complexifier, factorization and deformation, for construction of 
coherent states are presented for one dimensional nonlinear harmonic oscillator (ID 
NLHO). Since by exploring the Jacobi polynomials P,i'^'s, bridging the difference be- 
tween them is possible, we give here also the exact solution of Schrodinger equation of 
ID NLHO in terms of Jacobi polynomials. 

1 Introduction 

For description of the dynamics of a non-linear harmonic oscillator (NLHO) [T], we will 
encounter with three types of annihilation operators, their eigenfunctions can be considered 
as coherent states: (i) The annihilation operator, which is obtained from complexifier, 
introduced by Thiemann [2]; (ii) The annihilation operator, which is constructed out of 
the annihilation operator of simple harmonic oscillator (SHO), by using the method known 
in the nonlinear coherent states theory [3l 15], and (iii) The annihilation operator which 
appears in the factorization of the shape invariant potentials |5j. 

Our work is divided in two parts: The present contribution dealing with the mathe- 
matical construction of three types of coherent states of a nonlinear harmonic oscillator. 

In the second part, the relation between them, their properties, and some physical 
realization, will be given by exploring the properties of Jacobi polynomials (instead of 
Hermitian polynomials in the SHO theory). 

The paper is organized as follows: In section [2l by introducing a canonical transfor- 
mation, we give a prescription for quantization of classical Hamiltonian, by which the 
Schrodinger equation takes a familiar form of a particle in a tanh^-potential. Section [3] is 
devoted to solving the Schrodinger equation, and finding the eigenvalues and explicit form 
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of eigenfunctions of the Hamiltonian. In section U classical and quantum complexifier for 
this system are derived. Finally in sectional three types of coherent states related to the 
different annihilation operators will be presented. 



2 Schrodinger equation of ID NLHO 

The equation of one dimensional nonlinear harmonic oscillator is given by [T], 

m[(l + Xx'^)x — Xxx^ + uP'x] = 0, 



where uj is the angular frequency, A > is a real parameter with physical dimension 
{length)'"^ and for A = the equation of motion of an oscillating particle with mass m is 
obtained. The equation of the ID NLHO can be obtained from the following Lagrangian 
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(3) 



2 V 1 + Aa;2 
The canonical momentum is given by 

dL mx 

P 



dx 1 + Ax^ 



and from this equation we get 

X = {1 + Xx'^)p/m, p 



Xxp^ 



rauP'x 



m (1 + Xx'^Y 
Moreover, the Hamiltonian can be obtained as follows 



. (1 + AxV^ 

H = px - L = h 

2m 



2 2 

mui X 



2(1 + Ax2) 



(4) 



(5) 



Note also that p = —dH/dx. From the energy conservation H = E, H describes an 
oscillator with an amplitude dependent frequency. 



x{t) = Asm{nt + ^), A = 
Thus, by using equation ([T]), we get 



2E 



moj'^ - 2EX ' 



2E 
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w Vl + A^2 ' 
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~2X 
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rriLo'^ - 2EX 



m 



1 + XA^ 



Now, by define the following canonical transformation 

y/TTXx^p = P, sinh" ^ (VXx) = X 

V A 



(6) 



(7) 



(8) 
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we get 



X 



P 



m A cosh^VAX dX 2A 

And for these transformations, it is easy to check the following Poisson brackets: 



{x,p](x,P) = l = {X,P} 



{x,p) 1 



where 



lim X = X, lim P = p, 



dX 



ox op op ox 



dx 



The Hamiltonian in the new coordinates is given by 



H{X,P) = ^ + ^tanh'{VXX). 

ZTTi LA 



(9) 

(10) 
(11) 

(12) 



The form of potential energy in this coordinate system has clear relation to the force 
equation in ([9]), and also P = —dH/dX. 

Remark. Even though H looks like the Hamiltonian of a particle in a potential given 
by a tanh^ function (or tan-^-function if A < 0), but under the canonical transformation 
X — >• X and p — )• P, the relation between the kinetic and potential energy is more subtle 
and the former is not given simply by P'^ /2m. 

Before quantization of the Hamiltonian and writing the Schrodinger equation, we have 
to consider the following relations in Schrodinger representation: x — )• x, p — t- p = —ihd/dx 
and [x,p] = ih. Thus there are many representations for P regarding the Taylor expansion 



of \/l + Ax^. Two of these representations are given by, 

A = -iWi + Ax2-^, with pf = -r? ( \X^ + {I + \x^)-^\ 

ax \ ax ax^ ) 



and 



d 



-ih—y/l + Xx^, with P| 
dx 



{ X + 3Ax-^ + (1 + Ax 
* dx 



dx^ 



(13) 



(14) 



However the commutation relations of these two representations with X has the same 
result. 



[X,Pi]= sinh-^(VAx), -ihy^l + Ax2-^ 

V A dx 

[X,P2]= -^sinh-^(VAx),-«;i-^\/l + Ax2 

V A dx 



ih. 



Thus both representations can be realized as —ih[d/dX). But we should choose one of 
them from the outset, because these two representations do not commute and we would 
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have the ordering problem of operators as is usual in quantization theories. On the other 
hand the similarity transformation 

e/Wp^e-/W = p^^ f^x) = ln(l + Xx^)-^/^ (15) 

relates Pi and P2. And the Hamiltonians Hj = i^ + V{X) with j = 1,2 are also related 
by the same similarity transformation, 

e^(^')i?ie~-^(^') = H2. (16) 

Their eigenvalues are identical and their eigenfunctions are related by e-^^^^^pi = ip2- By 
using the first representation Pi = P = —ih{d/dX), we get the Schrodinger equation 

H^{X) = [-^^ + ^tanh\^X)MX) = E^{X). (17) 

3 Solution of the Schrodinger equation of NLHO 

The Schrodinger equation for tanh^-potential ()17p can be solved analytically by using 
the methods developed in [6l [7], and in this section we review the the solution for later 
references. 

3.1 Hypergeometric differential equation 

First by defining three dimensionless parameters 

r- rr?u? 2mE 

VXX = z, j^ = v, -^ = e. (18) 

Then the equation (jl7p can be rewritten as 

^-^^^ - [wtanh^z- e]v3(z) = 0. (19) 

Now, if we consider the following transformation 

1 1 — 2u 

z = tanh-^ -^=^= = sinh"^ i{l - 2u), (20) 

2 y/u{l - u) 

then we obtain u = (1/2)(1 + isinhz), and cfcp/dz"^ = -n(l - u){d^>^/dv?) - (1/2)(1 - 
2u){dip/du). Therefore we have the following second order differential equation, 

n2(l - uf^ + \u{l - u){l - 2u)^ - - 2uf + eu{l - u)] if{u) = 0. (21) 
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Next, we define a new function 

$ = n^/4(i _ ufl^if, (22) 
tlie first order differentiation will be removed from the equation and we get 



9. ^'yd?^ 



. 1 , 2 / 1 N 2) V 

{e + --v)u -{e + --v)u + --- 



^{u) = (23) 

which could be transformed to a hypergeometric equation by defining the following function 

y{u) =u-''{l-uy^, (24) 



where = \ + \\J\ + v and a"^ — a = (—3/16) + {v/'i). After some calculus, we get 

^(1 - + 2^(1 - 2^)^^ - - ' - 8 " 2^^^''^ = ^^^^ 

This has the desired form of a hypergeometric differential equation with well known solu- 
tions. 

3.2 The general solutions and energy eigenvalues 

The equation ()25p has the general form 

u{l-u)^ + [-/-ia + l3 + l)u]^-X/3w = (26) 

which is satisfied by Gauss hypergeometric functions for \u\ < 1. This equation has exactly 
three singular points: u = 0, n = 1, and u = oo. When none of 7, 7 — — /?, or a — /3 is equal 
to an integer, two fundamental solutions are known for u near each of the three singular 
points of the differential equation. Because of the Schrodinger boundary condition, we 
study the solutions near u = ±00. Two linearly independent solutions for \u\ — )• 00 is given 

by [8J, 

u-"2Fi{a,l + a-TA + a- u"'' 2i^i(/3, 1 + /3 - 7; 1 + - a; n"^). (27) 

There is an equivalent representation for each of the above solutions, which we give here 
for later applications 

n''-^(u-l)^-"-'^2i^i(l-/3,7-/3; a-/3+l; n"^), u'^-^{u-l)^-''-^2Fi{l-a,-f-a; /3-a+l; u'^] 

(28) 

The general form of 2-^1 has the following expansion 

J7 f h \ (a)n(b)n ^" 
2F1 (a, 5, c;n) = 2^—- — , (29) 

n=0 
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where e.g., (a)„ = a(a + l)(a + 2) • • • (a + n - 1) = [(a + n - l)!]/[(a - 1)!], (a)o = 1. The 
Gauss serie reduces to a polynomial of degree n in u, when a or /? is equal to — n, n S No- 
The differential equation (j25p is like ()26p with the following parameters: 



7 = 2c7, a + /3 + 1 = 4c7, a/3 = 2(T^ + e 



1 V 
8~ 2' 



Prom the last two relations one can easily find 

1 



a = 2a- ^- \/^eT^, (3 = 2a - - + y/-e + v. 



(30) 



(31) 



Because of the second relation in (pOj) , the sign of \Jv — t in expressions of a and /3 should 
be opposite, and it seems that they could be chosen conversely. But since some restriction 
will be inserted on /3, we have to choose the sign as given above. 
According to the relations (fM|) . 

(32) 



<^(z)=n-l/4(l_^)^-l/4)y(^)^ 



which is defined as before by 

n = -(1 + i sinh z) = -(1 + ig), n(l — n) 



\u\ 



(l + g2)/4= [cogh2z]/4 (33) 



so that 



<^(^) = n'^-V4(i_^)<x-l/4) 



Fi (a, 1 + a — 7;1 + Q — /3;ii 



+02-^-^ 2i^i (/3, 1 + /3 - 7; 1 + /3 



a; u 



If 2; — )■ ±00, then g — )• 00, w ^ — )• 0, and we have 2Fi{u ""^ = 0) = 1, thus 



(34) 



lim (p{z) 
2— > it 00 



u 



,(T-l/4) 



(35) 



If we absorb temporary the phase factor of n~", in ci, C2 and set \u\ instead of n, then 
the boundary condition will be satisfied for u = —00 and u = +00 simultaneously. Now, 
by using the definition of a in ()24p and its relation with a, /3 in (j3ip we get 



1, 



liin ^(^) = ci^[yiT?]-(^) + c2^[v'r+7](^). 



(36) 



The second term goes to infinity as z — >• ±00, so we have to set C2 = 0. The serie in ([3^ 
converges if < 1, or equivalently, if(l + g^)/4>l, then \q\ > ^/S. The requirement 
for the wave function being everywhere finite, forces us to the analytical continuation of 
(p{u) in the interval \u\ < 1 or (7 in (— -y/S, \/3), (the endpoints are not included). The idea 
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is to define a transformation to map as large a region of the complex plane as possible onto 
the interval ^(1 - i\/3) < n < i(l + iy/3) or (e'^'^/^ < u < e*'"/^). If arg(-?i;) < vr, then 
the analytic continuation of 2-^1 (^^) is given by |9j 



iFi{a,b; c;w) 



vr 



where 



T(c)(-w)^°- 1 
■ ( ^^ rfh\rt . 2Fi{a,a - c+ l;a - b + 1; -) 

sm(7r(o — a)) Li [b)L [c — a) w 

rffj^^'K 2F^{b,b-c+l;b-a + l--] 
r(a)r(c — b) w 



r(7?)r(i-r?). 



(37) 



sin(7r7/) 

After setting C2 = in and a = a, 6 = a — 7 + I, c = 1 + a — /3 and = in (j37|) . 
we get 

ifiz) = (-l)°ciu"-i/^(l-'u)"-i/'*) 2i^i(a, /3; 7; ^x)-^^^-^ 2i^i(a-7+l, /3-7+1; 2-7; , 

(38) 

where 

A = r(l + a-/3) 

sin 27ra r(a - 7 + l)r(l - (3) ' ^ ^ 

^ r(l + a-/3) _ ^^^^ 
sin 27r(T r(a)r(7 — ^) 

We have to handle two problems in the above transformation: First, due to the presence 
of r(l + a — /?) in the numerators and sin27rcr in denominators of A and B, the cases 7 = 
2(T G Z and a — /? G Z cannot be handled using these formula. Second, this transformation 
cannot be used for u = exp(ibi7r/3) = (1/2)(1 it i^/S), because these points are mapped to 
themselves or each other, and at these points the serie diverges. 

By considering (I24p and (I30p . the first problem means that v has to be an integer and 
this is not the case for general physical parameters. In such an special cases, there is 
other transformation to tackle the problem [9j. But the second one is the major problem. 
Since the wave function is finite at these two points, the serie has to be terminated. The 
first possibility is that a or 6 in the serie representation of 2F1 in (j29jl would be equal to a 
negative integer. In the denominator of A and B, if the functions r(t)'s have a pole position 
at one of its arguments, then the exponentially growing term is absent. If a = a = — n 
then B = 0, since it has the T{a) in its denominator. In other hand if we set (3 = n' or 
1 — 13 = —n, then A = for the same reason. Both possibilities have the same result as we 
will see below, i.e, if one of A or i? vanishes, the other serie terminates. Now, if we choose 

1 - /3 = -n, n G No. (41) 
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then, ^ = and by using T{z + n) = (z)„r(z), T{z — n){z — n)n = T{z) and T{z)T{l — z) = 
vr/ sin(7r2;), we get 

r(-2. + 2)(-2. + 2) (42) 

^ ^ (4cj-2-2n)„ ^ ^ 

From and the condition /3 = n + 1, we can find the discrete energy eigenvalues, 

2 



- + v - (n+ -) 



(43) 



so that there is an upper bound for n, 



n + -< Vv - er 



(44) 



(This condition is satisfied for finite values of A and for the special case A — )• oo, t; — )• 
and we have = —n?). Since for an atomic system v is relative large number, the energy 
can be considered always positive. If we use the physical parameters given in (llSp . then 
we obtain 



^" 2A 2m 
And in the limit A — )• 0, En becomes 



1 m?uj'^ 1 
4 + W-^"+2) 



(45) 



lim En = (n-\ — )tiuj 
A^o 2 



(46) 



as we have expected. 



3.3 Eigenfunctions of the Hamiltonian of NLHO 

Now we consider the equivalent representation of the solutions given in (j28p for the differ- 
ential equation ([25]) with parameters given in (f30|) and C2 = 0, 



ip{z) = ciu''-^/\l - uf-'^/^^u'^-^{u - ly-^-l^ Lfi{1 - /3, 7 - /3; a - /3 + 1; -u-i)] . (47) 
By using the same transformation as given in (|37p . we get, 



2F1 (1 - /3, 7 - /3; a - /3 + 1; u-') = T{j - l)r(2 - 7)r(a - /3 + 1) 



Lr(7-/3)r(a)- 



;Fi(l-/3,l-a;2-7;n) 



(_^,-i)-{7-") 
r(l-/3)r(a-7 + l) 



2^1(7 - /3,7 - a;-/;u) 



The condition 1 — (3 = —n, implies r(l — /3) — t- 00 and the second term in the bracket 
vanishes and the first hyper geometric function reduces to a polynomial of degree n. The 
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coefficient of the first term in RHS is again B, as given in ()39p . so that we have the following 
eigenfunctions belonging to the eigenvalues En- 

^^{z) = -(-l)"ci5u'^-i/^(l - uy-^/^\l - 'u)T-"-^u^-^2i^i(l - /?, 1 - a; 2 - 7; u). (48) 

The polynomial representation of 2-^1 is given by [8], 

2Fi{-n,a + 1 + b + n,a + l;u) = Pi°'^)(l - 2-u), (49) 

(Cl + i)n 

where plf'^^ is Jacobi's polynomial, 

^ ' m=0 \ / \ I 

By comparing to the relations ([l8|) . ([30]) . ([32]) and the condition 1 — /3 = — n, it is easy to 
see that: 1 — 2n = i sinh z, a=l — 7=1 — 2cr, 6 = 7— (a; + /3) = l — 2cj, and by using 
the relations (pOj) . (jl3|) . we have a = 4a" — n — 2. Therefore the eigenfunctions (j^Hjl can be 
written as 

(^„(z) = (-l)4--"-icii?n--+t(l - n)-+l p(-^-+^'-^-+^(-isinhz). (51) 

(-2fT + 2j„ 

Thus in the old coordinates the eigenfunctions are given by 

^.(x) = (-l)^— ^cii3(l + Ax^)-+t p(-^-+^'-^-+^)(-.VAx). (52) 

(-2(7 + 2j„ 

and by Rodrigues' formula for special case \f\x G M 

In terms of physical parameters, B is known and we could choose c\ in such way that the 
eigenfunctions become normalized. Now we calculate X\m.\^QLpn(x) in several steps: 

• From (p^ we have lim;\_j.o o = 2X62' ~- 2 ^ ^' '^here }? := h/moj. So that lim;^_j.o(l ~ 
2(t) = — ^. From now on we use the limA_s>o ^-nd lim^^oo alternatively. 

• We note that u'^^iil - uy+i = (l/4)-'^+t [1 + sinh^ \/\XY''+'i. For large values 
of ^, we have sinh^ ^/XX ~ Ax^, (thanks to the canonical transformation between X 
and X as given in (jSj), so in this limit, 

(l + Ax^)-«/^ = (l + ^^)-^/2 = e-^/^^^ 
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There is a relation between Jacobi and Hermite Polynomials [TO] 

where C ^ &iid we have the same result for ^ — )• — oo. Therefore if we set C = 
and considering that, 



then for large values of ^, we have 

on J, I ™ 

lim — = F„(x/6). 

5^00 (-e)("/2) " 

In the limit A — t- 0, by choosing the appropriate form for ci, we obtain the normalized 
eigenfunctions of the Hamiltonian of simple harmonic oscillator, 

^n(x) = , e-^'/^'^HrA (53) 

4 Complexifier method 

In this section we define complex coordinates on phase space and their corresponding 
operators by using Thiemann's method [2]. 

4.1 Classical complexifier 

The classical complexifier is defined to be kinetic energy function divided by a constant of 
dimension frequence w, which can be expressed as 

kinetic energy 

complexiner = . (54 j 

OJ 

We construct complex- valued function "z", corresponding to the position coordinate x, by 
taking the position functions x and applying repeated Poisson brackets with the complex- 
ifier, 

' n=0 

where Poisson bracket is given in pTj) and {x,C}(o) = x, {x,C}(„_|_i) = {{x, (7}/„), C}. 
According to the original classical Hamiltonian ([5]), the kinetic energy is m) 
so that the complexifier is defined by 

C=^l±^. (56) 
2moj 
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Some of the leading terms of Poisson bracket are given by, 



so that 



2mu 



[2p{l + Xx% {x,C}(2 



{x,C} 



^2n 



(2n) 



(2) 



:[Va(1 + Ax2)x], 



[muj 



\2n 



{x, C}(2„+1) 

Then (15511 can be calculated, 



^ -A"(l + Ax2f+i. 



(ma;)2"+i' 



cos 



iy/{l + Xx^) . /pVA(1 + Ax2) 

X A — = sm — 



muj 



Vx 



In the limit A — )• we have 



lim z 



X H — p. 

2 y2mijj 



(57) 



(58) 



The physical dimension of z, z* is {angular momentum)^^^'^\ so that their Poisson bracket 
is dimensionless. They are known as complex coordinates on phase space and we also have 



limujz*z=p /2m + muj x /2. 
A— >-o 



(59) 



It is important to note that when A 7^ 0, this transformation is not canonical, i. e., in spite 
of {x,p} = 1, 

'py/X{l + Xx^y 



{z,z*} 



Xx^ + cos 



(60) 



and limA-!>o{-Z) -2*} = —i- 

We can find the corresponding complex coordinates for the canonical transformed X, P, 
according to dS]). In the theory of canonical transformation in classical mechanics we have 
the following statement: 

If (p : — )■ M^-'^ is a canonical transformation, then, 

{h o (j), g o (f)} = {h,g} o(j), 

where / is degree of freedom, in our case / = 1. Accordingly, 

1 P 

{x 0, Co = {x,C}(„) o0, (l){x,p) = (— sinh\/AX, — ; — 7=—). 

VA coshyAA 
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Our new complex coordinates is given by 
I °° f\ 



oo 

' n=0 




2X 



sm cos + i cosh(VAX) sm 

\ muj I \ mijj I 



(61) 



Therefore we get to the complex coordinate related to x, 



Z = J^sinhiVXX+'-^], Z* = J'-l^sinhiVXX-'-^]. (62) 
V 2A \ muj I V \ J 

In a similar way we get to the complex coordinate related to X, by using the same com- 
plexifier, 

oo 



n=0 



and we also have 



A=J y + , A* = J X ; . (63) 

V 2 V 2 ^2^ ^ ^ 

It is seen that Z and A have a similar relation, as x and X, 



and there is a similar relation for A* and Z* . 

Remark. In comparison, these relations are counterparts of the coordinates x and X 
as given in ([8]). Therefore Z is complexification of x, such as A is complexification of X, so 
that one can say, the complex coordinates are related according to the following diagram, 

X Complexifier Z 

ct I I 
X Complexifier A 

where ct stand for canonical transformation. By referring to ([9]) we can see, 



A = {A,H} 

2 ym cosh^ \/AX y \^ ' ' muj 

{Z,H}. (65) 



• Imuj I P iuj s\VL\\^/\X \ , ( rr^r iV^P 
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The Poisson bracket of Z, Z* leads to the same result as (|60|) , 



{Z,Z*}(^x,P) - IT 



cosh(2V AX) + cosh(- 



muj 



Xx^ + cos^ 



pVA(1 + Ax2) 



-i cosh 



2A 



A cosh 



muj 



2A 



A* 



— {Z,Z*](^^^p) — -i{Z,Z*](^AA*) 
{A,A*}^x,p) = -i = {vl,^*}(x,p)• 
4.2 Quantum Mechanical Complexifier 



(66) 



Analogously, the complexifier in quantum mechanics is defined by kinetic energy operator, 
so that the corresponding operator of complex coordinate is given by 




2 ^ n\ (ihY 

n=0 



(n) _ ^-C/h^^C/h 



(67) 



where z plays the role of annihilation operator and as given in ()13p . 



C 



2muj 



1 . d 



+ :tAx 



2 dx 



(68) 



and [x,C'](o) = x, [x,C](^n+i) = [[^) C'](n); C*]- The calculation of different order of com- 
mutator in (x,p)-coordinate system is cumbersome. It is easier to go to the canonical 
transformed coordinates X, P. But again we have to be careful in definition of kinetic 
energy operator. 

As it is done in the case of classical complexifier, the equivalent version of (I67p is given 

by 



Z 



OO , .^ 




2 ^ (i/i)"n! 

n=0 ^ ^ 



—= sinh \/A-'^ , o 

^ 2ma;dX2 



(n) 



n=0 



r^N n (_1)« 



n! 



sinh \/AX 



72 n 



(n) 



(69) 



where rj = — . The following commutation relations will be useful in sequel, 



sinh-v/AX , 
cosh VXX, 



dX^ 



-A sinh VXX - 2\/A cosh VXX 
-A cosh VXX - 2\/A sinh VXX 



_d_ 

dX' 

_d_ 
dX' 
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Therefore we can calculate the repeated 
terms in summation read as 

Sn = 



commutators more easier. Some of the leading 

So = sinh VXX 

VXX)X + (cosh VXX)2VX — 
dX 



Si = (^) ^(sinh\/Xx)A + (cOBhVAX)2VA^^ 
2 2! 



S3 = C^)'':^ (^(sinh^X)(A=' 



2 V dxj 

(-f— ( (sinh VAX)(A^ + 4A-^) + (cosh VXX)(4A^/^ — ) 
2 2! \ dX^ dX 

ihV\X)(\^ + 12A^ ^) + (coshV\X)(6x'^^'^ h 8A^' 



+ (cosh^/AX)(6A''/^ — 



h SA-"'^ : 

dX dX' 



dX^' 



i(A'' + 



24A'' + 16A 

dx2 ^ 



^-^) + (cosh VAX)(8A''/^— +32A^/^-^-^) 
dX* dX dx3 _ 

(-)■''— f'(sinh^X)(A^ +40A*'-^^— + 80A-''-i— ) + (cosh ^/AX)(10A^''^— +80A^''^-i^ 
2 5! \ dX2 dX<l dX dX-' 

2 6! y 



(^)^i (^(BinhVAX)(^^ 

S5 = (^r^^((.i 



d^ „,„ d^ 

h 32A^' ^ I 1 

dX3 dX5 7 

H eOA*^ + 240A* +64A^ -) + (cosh V^X)(12A^^^^ h 160A^^^ + 192A''^^ ; 

2 61 V dX^ dX* dX<i'^ dX dx3 dXf 

S7 = {-)''— ('(sinh VAX)(a'' + 84A''-^ + 660A^-^ + 448A''-^^) + (cosh VAX)(14A^^/^ — + 280A"''^-!i— + 672A''/^-^ + 128A^/^-^^) 
^ ^2' 7! V dX^ dX* dXS'^^ dX dX3 dX^ dXT' 

Then after summation over appropriate terms we get, 
Z = 



) 



,r,A/2 



2A 



sinh(\/AX) cosh(r/\/A--^) + cosh(\/AX) smh{'q\f\-^] 
dX dX ' 



2A 



sinh(\/AX) cos( ) + i cosh(\/AX) sin( ) 

moj moj 



(70) 



Here we have to note that again, like the classical complexifier, Z is complex operator 
corresponding to the x. If we construct complexification of X by using C we get. 



A = 



imoj 



X, 



X 



(n) 




2^ V2mujh' 



(71) 



where we have nniltiplicd the factor 1 / v^, to make A a dimensionless operator. By using 
the Baker Campbell Hausdorff (BCH) formula 

^{A+B) ^ ^A^B^-[A,B]/2^ B]] = = [B, [A, B]] 

and writing the hyperbolic and trigonometric functions in exponential representation, 



Xh/inm 

2A 



smh{V\X + 



mco 



2A 



sm 



^^/xx 



muj 



(72) 
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we can see the correspondence between classical, (|62]) . and quantum complex coordinates. 
Although in the quantum case, by considering the operators ordering problem, we were 
careful in using the trigonometric combination formula but thanks to the properties of 
these functions, we see this correspondence, up to a constant factor, nicely. The result 
commutator, then, will be achieved, 



2A mcj 



cosh(2\/Al) + cosh(-^^'^^ 



(73) 



Again, in the limit A — t- 0, 



hm [Z, ZT] = h, 

A— s>0 



and by referring to the equations ([60l [66l 



^mJ-[Z,P] = {Z,Z*} = {z,z*}, 
A-5>o in 

and 

i(ZZt + Z^Z) = ^e^^'" cosh Xb^ ( cosh^ VXX - cos\VXb^P/h)) . (74) 
2 2a V / 



According to the Hamiltonian (p!7[) it easy to see that 

A=^[A,H], ^ = hz,H]. (75) 
in in 

Again we see the correspondence between quantum and classical, (j65p . equations of mo- 
tions. 

5 Three types of coherent states 

5.1 Type 1: Complexifier coherent states 

At first, one can immediately recognize that by using (j7ip . 



hi 1 ^ 

— (iti + lit) = — + -mu^X^, and [A, it] = i. (76) 

2 '2ty\i 2 

It is easy to show that 

AMx') = 0, MX') = (^) ' «^p(-5)' C^^) 

in the limit A — )• we get the ground state of SHO. In the old coordinates 

M^)=(^ -P ^2 ' (78) 
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where 6^ = fi/{muj), and for every 7 G C, we can easily construct other eigenfunctions of 
A, by using the displacement operator, 

(X'|7) = e-l^l'/2e7^Vo(^')- 
Now, we have the standard coherent states in {X, P)-coordinates, 



'$57, 



(X)^ = \/2 6 K7, {P)j = V2 hmuj ' 
and it can be shown A\^) = 7I7). We have also the following relations, 



A =: = e^^' sinh[6V2Ai], A = sinh-i[e"^^'7(2AM^^]- 

which are similar to (j64p . We have divided both sides by Vh, for dimensional convention. 
Then 

Z'h) = /^e"'' sinh[6x/2A7]|7> = /(7)|7> 

To find the explicit form of the eigenfunctions of Z, we assume that \z) is the desired 
eigenket, Z'\z) = z\z). Using the (over) completeness of |7)'s, Tr~^ f c?^7|7)(7| = 1, and 
the inner product of standard coherent states |(^;|7)P = exp{—\z — 7p), z is given by the 
following equation, 

z = TT-^ / (iV(7)K^l7)P = / d27sinh[6V2A7]e-l^-^l 



5.2 Type 2: f-deformed coherent states 

Now we introduce another form for energy eigenvalues (I45p . We have to find a n-dependent 
function /(n), so that En is written as, 

i5;„ = ^[(n + l)|/(n + l)p + n|/(n)|2]. (79) 
This function is given by 
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According to the upper bound condition (jM]), we can consider the /(n) as a real function; 
and we have 

lim f{n) = 1, hm En = (n + -)nuj. (81) 

A— >0 /(n)— >1 2 

First we give the following representation for the Hamiltonian in Fock space of SHO, by 
using the relation (|l5]) . 



H = ^(bb^ + b^b) 

= ^[(n + l)|/(n + l)p + n|/(n)|2] 



- + v-{n + -) 



(82) 



Then, we can formally construct the relations between creation and annihilation operators 
of NLHO and SHO, 

b = af{n), b^ = f{n)a), a) a = n, (83) 

so that 

T 



[b, 6t] 



[(n + l)|/(n + l)|2-n|/(n)p 



- +t; - (n+ -) 



(84) 



The coherent states are given as right hand eigenfunctions of the new annihilation operator 

b\f,l3) = l3\f,/3). (85) 



5.3 Type 3: Generalized displacement operator 

The Hamiltonian (I17p can be easily factorized in the following way. 

1 



(86) 



where 
B = 

and 



+ 



2A 



tanh(\/AX) 



-h d I mu!^ 
+ 



/2m 



2A 



tanh(VAX) 



' ^ cosh2(VAX) 1 + Ax2- 



(87) 
(88) 
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The ground state equation for annihilation operator BiPq{X) = 0, leads to 



iPo = ci[cosh(^/AX)](-i/^''') = ci[Vl + Ax2](-i/A6^). 
If we define ^ = 1/(2A6^), so that for A — )• 0, — )• oo, then we have 



(89) 



1 x'^ 

lim ipoix) = lim ci(l + = cie 

A-s>0 5-!>oo ^ 20^ 



This is the ground state of simple harmonic oscillator. Since this system belongs to the 
shape invariant class, the ground state wave function should be of the form 4>o{X) = 
C2[cosh(\/AX)](-"-^o), so that 



HMX) = EoMX), Eo 



2m 



'2 ^ V 4 ^ A264 



with a 



2m 



where (t>o{X) and Eq are exactly the same as we have in ([52]) and (|15|) . Again, for A — )■ 0, 
£"0 — )• {huj/2), and \(po) — )• |0), i.e., the ground state of SHO. Now, we construct the 
generalized displacement operator. 



And the coherent states are given by 

D{CMo) = \0, and ^10 = CIO- 



(90) 



(91) 



Conclusion and final remarks 

Three types of annihilation operators and their eigenfunctions, i.e., corresponding coherent 
states, are introduced for ID NLHO. All of them are dependent on a deformation parameter 
A, which can be realized as a control parameter in different physical problem (as we will 
see in the second part of this work), and for A — )• 0, they are identical with the annihilation 
operator of SHO. 

The most important problems, we are going to address in the second part of this 
work, are: Presentation of the explicit form of coherent states and investigation their 
most important properties, such as resolution of the identity, continuity in label (since 
the coherent states will be labelled by points of the phase space), and the relationship 
between these coherent states, by exploring the properties of Jacobi polynomials P"'^ as 
eigenfunctions of the Hamiltonian of NLHO. They play the role of Hermitian polynomials 
Hn in the theory of SHO, so that the corresponding coherent states may be constructed by 
their superposition, in a similar way as the standard coherent states of SHO are constructed 
out of HnS. 
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